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ABSTRACT 

New boron abundances for seven main-sequence B-type stars are determined 
from HST STIS spectroscopy around the B III 2066 A line. Boron abundances 
provide a unique and critical test of stellar evolution models that include ro- 
tational mixing since boron is destroyed in the surface layers of stars through 
shallow mixing long before other elements are mixed from the stellar interior 
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through deep mixing. The stars in this study are all on or near the main-sequence 
and are members of young Galactic clusters. They show no evidence of mixing 
with gas from H-burning layers from their CNO abundances. Boron abundances 
range from 12+log(B/H)<1.0 to 2.2. The boron abundances are compared to 
the published values of the stellar nitrogen abundances (all have 12+log(N/H) 
<7.8) and to their host cluster ages (4 to 16 Myr) to investigate the predictions 
from models of massive star evolution with rotational mixing effects. We find 
that the variations in boron and nitrogen are generally within the range of the 
predictions from the stellar evolution models with rotation (where predictions for 
models with rotation rates from to 450 km s _1 and //-barriers are examined), 
especially given their age and mass ranges. 

Three stars (out of 34 B-type stars with detailed boron abundance determina- 
tions), deviate from the model predictions, including HD 36591, HD 205021, and 
HD 30836. The first two of these stars have much larger boron depletions than 
are predicted for their spectroscopic masses and very young ages, even adopting 
the highest rotation rates from the model predictions. HD 36591 also shows no 
significant nitrogen enhancement, as uniquely predicted by the rotating stellar 
evolution models. HD 205021, however, has a small nitrogen enrichment which 
could be explained by stellar rotation or mass transfer since it is in a binary 
system. The spectroscopic mass for the third star, HD 30836, is marginally lower 
than expected given the rotating model predictions for its age and boron abun- 
dance. This star also has no significant nitrogen enhancement, thus even though 
it is in a binary system it does not show the nitrogen enrichment expected if it 
has undergone mass transfer. Therefore, the results from these three stars sug- 
gest that rotational mixing could be more efficient than currently modelled at 
the highest rotation rates. 

Subject headings: stars: abundances - stars: evolution - stars: rotation 

1. Introduction 

Rotation is recognized as an important physical component in understanding the evo- 
lution of massive stars and yet is a theoretically challenging problem. Rotation affects the 
lifetimes, chemical yields, stellar evolution tracks, and the properties of supernova and com- 
pact remnants (Heger & Langer 2000; Maeder & Meynet 2000, 2005). The new rotating 
stellar evolution models also address long standing problems such as the origin of the B[e] 
and WNL/Ofpe (slash) stars, the distribution of red to blue supergiants on the HR dia- 
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gram, and the unseen post main-sequence gap predicted in all standard stellar evolution 
scenarios. With rotation, it is also possible to explain the variations in the surface helium, 
carbon, nitrogen, and oxygen (He and CNO) abundances in OB stars on and near the main 
sequence (e.g., Gies & Lambert 1992, Herrero et al. 1992, Cunha & Lambert 1994, Den- 
nisenkov 1994, Lyubimkov 1996, Lyubimkov et al. 2004). Several of these observations have 
been used to constrain the various mixing prescriptions used in the new models of massive 
star evolution with rotation. Additional observations are now necessary to test the model 
predictions and provide new constraints for the transport of angular momentum and chem- 
ical species in rotating massive stars. One avenue of observational testing is to determine 
the helium and CNO abundances in massive stars in young star clusters in the Galaxy and 
Magellanic Clouds (Evans et al. 2005). This tests the metallicity effects predicted by the 
models. Another line of research includes the determination of the light element abundances 
at the surface of massive stars. This tests the earliest stages of rotational mixing (timescales 
and mixing efficiencies), before hot gas from the stellar interior is observable at the surface 
(Fliegner et al. 1996, Venn et al. 2002). 

The abundances of the light elements, lithium, beryllium, and boron (LiBeB) are known 
to be sensitive to rotational mixing in both low and high mass stars. LiBeB is destroyed on 
exposure to protons at temperatures too low for H-burning to have occurred (<6xl0 6 K), 
therefore even shallow mixing at the stellar surface induced by rotation can lead to LiBeB 
depletions. In low mass stars, variations in the surface abundances of all three elements 
have been traced to rotational mixing (e.g., Boesgaard & Lavery 1986; Pinsonneault 1997). 
In high mass stars, only boron 1 has been available. Spectroscopy of the B III feature at 
A2066 using the International Ultraviolet Explorer (IUE) archived spectra or the Hubble 
Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) or Goddard High 
Resolution Spectrograph (GHRS) have made boron abundance determinations possible in 
B-type stars. Significant variations in the boron abundances in hot, massive stars have been 
observed (e.g., Proffitt & Quigley 2001, Venn et al. 2002). 

Boron depletions in B-type stars are predicted to be associated with nitrogen enhance- 
ments. This is because mechanisms that can deplete boron at the surface of a star will also 
mix the surface with CN-cycled gas, i.e., gas from H-burning layers where the CN-cycle has 
converted carbon into nitrogen. This is true whether boron is destroyed through rotational 
mixing (where hot CN-cycled gas from the interior is mixed to the surface) or through bi- 



1 There are no readily accessible spectral lines of Li in hot stars, and the resonance lines of Be II near 3130A 
are predicted to be very weak in B-type stars, and they occur at the atmospheric cutoff, thus these spectra 
would be difficult to access from ground based telescopes. Currently, there are no published observations of 
the Be II lines in B-stars. 
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nary mass transfer (where CN-cycled gas is deposited on the surface of the star; Wellstein 
2001). Simple initial abundance variations and mass loss from B-type stars are ruled out 
(see Venn et al. 2002). However, there is one unique characteristic of the rotational mixing 
scenario that makes for an interesting and exciting new observational constraint. Boron 
depletion occurs before nitrogen enhancement. This is because rotational mixing taps the 
surface boron- free layers first, then subsequent deeper mixing can penetrate layers where 
H-burning has occured via the CN-cycle, converting carbon into nitrogen. Thus, the initial 
phases of rotational mixing in stars is revealed by depletions of boron, followed only later by 
nitrogen-enrichment and carbon-depletion, a different signature from binary mass transfer. 

Some B-type stars with boron depletions, but no nitrogen enrichments, were found by 
Venn et al. (2002) and intrepreted as concrete and unambiguous evidence for rotational 
mixing on the main sequence. However, two stars with masses of 12-13 M Q showed unchar- 
acteristically large boron depletions. These two stars required models with higher masses 
(20 M ) and the highest rotation rates (~450 km s" 1 ) to reproduce the boron depletions 
(from Heger & Langer 2000 models). Since very few stars are expected to rotate at these 
high speeds, and because of the mass discrepancy, this suggested that rotational mixing may 
be even more efficient than currently predicted. 

Presently, 32 solar neighborhood (see Table 9 in Venn et al. 2002) and two stars in 
the Small Magellanic Cloud (Brooks et al. 2002) B-type stars have boron abundance de- 
terminations. Of these, only nine have been determined from high quality HST STIS or 
GHRS spectra of the relatively unblended B III 2066 line. In this paper, we present new 
boron abundances for seven additional B-type stars from HST STIS spectroscopy. These 
stars were selected from the group of B-type stars examined with International Ultraviolet 
Explorer (IUE) that may have low boron abundances (Proffitt & Quigley 2001). They are 
also all in young galactic clusters (for age estimates), and they show no nitrogen enrichments 
(12+log(N/H) <7.8). Therefore, the determination of boron abundances in these stars pro- 
vides a unique observational test for studying the earliest phases of massive star rotational 
mixing effects. 



2. Target Selection and Observations 

Seven B-type stars were observed using HST STIS spectroscopy of the B III feature at 
A2066. Previous analyses of the target stars (Gies & Lambert 1992; Cunha & Lambert 1994) 
show that the targets have normal, unenriched nitrogen abundances in their atmospheres. 
Analyses of IUE spectra by Proffitt & Quigley (2001) indicate that these targets have very 
weak or absent B III lines, however the resolution of IUE and signal-to-noise of some of 
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these spectra make quantitative boron abundances highly uncertain. These 7 stars were also 
selected from young galactic OB associations in order that effects of rotationally induced 
mixing could be examined with age. Data considered in the selection of these targets is 
listed in Table 1. 

Spectra for six of the seven targets were gathered using the E230M grating and 0.2x0.05ND 
slit with a neutral density filter (ND=2) to avoid the MAMA-NUV brightness limits. This 
grating choice is sufficient to resolve the B III line in all of our targets since the instrumen- 
tal broadening of this grating is 9 km s _1 (given our slit choice), whereas the targets have 
rotational velocities <40 km s -1 . To attain the best possible signal-to-noise spectra, three 
different configurations of the E230M grating were used, each with a slightly different central 
wavelength (A c of 1978, 2124, and 2269 A) so that the B III line landed on different regions of 
the MAMA detector and subsequent combinations could better reduce pixel noise variations. 
Echelle observations with the E230M grating included 40 orders; only those orders in the de- 
sired wavelength region (AA2044-2145) were extracted and combined using standard IRAF 2 
packages. Spectra were then normalized using a low order spline3 polynomial; wavelengths 
listed here are in air and at rest. Exposure times, central wavelengths, and observation dates 
are listed in Table 2, as well as the signal-to- noise ratio of the final extracted spectra. 

The reductions for two stars, HD 36285 and HD 214993, were conducted using procedures 
different from those described above. HD 36285 (V=6.32) was faint enough to be observed 
with the higher resolution E230H grating (31x0.05NDA slit), and HD 214993 is a known 
variable star. The observations and reductions for these two stars are described below. 



2.1. HD 36285 

Observations of HD 36285 were done using the STIS E230H grating at the A2013 central 
wavelength setting. This star was too bright to be observed with any of the clear STIS 
apertures, but since we desired very high S/N in order to study the 10 B/ 11 B isotope ratio 
(Proffitt et al, in progress), we adopted the unsupported 31X0.05NDA slit, which has an 
attenuation factor of about 4. The global count rates were 1.2xl0 5 to 1.5xl0 5 cnts/s per 
observation, which is modestly under the 2xl0 5 cnts/s upper limit allowed for the NUV- 
MAMA detector. Using a long slit and dithering our subexposures in the cross-dispersion 
direction along the length of the aperture allowed for higher S/N. Nine different dither 



2 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation 
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positions with a spacing between positions of about 0.15" (or about 5 pixels) were used. 
This averages any fixed pattern noise of flat-fielding errors, and removes any bad pixels or 
detector defects. 

Since the 31X0.05NDA aperture is not a supported aperture, the standard reference 
files supplied by STScI cannot be used to produce extracted and flux calibrated ID spectra. 
While the measured throughput curve of the 31X0.05NDA aperture was used for calculating 
the scattered light background, we adopted the parameters supplied for the 0.3x0.05ND 
aperture. In addition, for some of the more extreme dither positions, the reference file which 
specifies the default locations of spectral orders on the detector had to be edited so that the 
calstis software could correctly identify each order. 

In principle, the use of a long slit with the echelle could increase the amount of interorder 
scattered light and cause an error in the background subtraction. There are some modest 
S/N (15 to 20) E230H observations of HD 39060 (f3 Pic) which were used to check for such 
effects (Program GO-7512, PI Lagrange; see Roberge et al. 2000). Observation o4g002020 
was taken at the 2013 central wavelength setting with the 0.1X0.03 aperture, while o4g002060 
was taken at 2263 with the 31X0.05NDA. Four echelle orders between 2125 and 2155 A are 
included in both observations. Uncertainties in the blaze function alignment at different 
central wavelength settings, and possible vignetting caused by the different locations of 
the orders on the detector in the two observations prevent an exact comparison, but when 
we reduce the 31X0.05NDA aperture observations of HD 39060 in the same way as for the 
HD 36285 data, these observations show good agreement for the extracted fluxes and suggest 
that any error in the background subtraction amounts to no more than 4% of the mean flux 
level. These observations also suggest that near this wavelength, the throughput of the 
31X0. 05DA slit is about 20% lower than its tabulated value. 



2.2. HD 214993 

HD 214993 (12 Lac) is a complex multi-mode (3 Cephei variable, with at least six distinct 
pulsational frequencies (Jerzykiewicz 1978) with periods between 0.095 and 0.236 days. The 
interaction of the different modes leads to considerable variation in both the amplitude of 
the pulsation and in observed profiles of absorbtion lines. Periods and amplitudes of these 
modes were measured by Mathias et al. (1994) using the Si III triplet at 4553, 4568, and 4575 
A. He fit the velocity curves with a sum of sinusoidal functions, although this fit does not 
perfectly track the velocity curve, possibly due to significant non-linearities in the pulsation 
amplitudes. 
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Our 11 subexposures of this star spanned about 0.195 d, or about 1 full cycle of the 
largest amplitude mode (0.193 d). We measured velocity shifts for each of the individual 
sub-exposures by cross-correlating the data in order 99 (~ 2045-2080 A) with our adopted 
synthetic spectrum for this star. The strong Zn II ISM line was excluded from this cross 
correlation. The measured velocities are plotted vs. observation time in Fig. 1. 

The 14 year span between the Mathias et al. observations and our STIS data is much 
too long to allow us to predict the phases of individual modes. However, for illustrative 
purposes we have overplotted a velocity curve that adopts Mathias et al. 's mode frequen- 
cies, amplitudes, and system velocity, but with phases for the individual modes arbitrarily 
adjusted to fit our observed data. The rms residuals of this fit are only about 2 km s -1 . 

In Fig. 2 we show for each subexposure the observed spectrum in echelle order 99. 
For our analysis we only used the data from subexposures 4, 5, and 6 which show the 
narrowest line profiles, which were taken when the star was near its maximum radius. We 
compare the coadded data for these subexposures to the coaddition of the other subexposures 
in Fig. 3. Note that the optical analysis of Gies & Lambert (1992) did not comment on 
the pulsational phase at which their observations were taken. This may introduce some 
additional uncertainties into their results for this star. 



3. Abundance Analyses 

Elemental abundances have been determined from LTE spectrum syntheses and AT- 
LAS9 model atmospheres (Kurucz 1979, 1989). Solar metallicity models were used through- 
out. The stellar temperature and gravity values were adopted from the literature; see Table 
3. For HD 36285, the temperature is from Cunha & Lambert (1994; CL94). For all other 
targets, temperatures are from Gies & Lambert (1992; GL92) scaled down by 3.4% to bring 
them into agreement with the CL94 temperature scale. 

LTE spectrum syntheses were constructed using the program LINFOR 3 . The spectral 
line list covers the region from 2044 to 2145 A, and was adopted from Brooks et al. (2002). 
This line list originated from the Kurucz (1993; CD-18) line list, including light elements, all 
lines in the iron-group, and heavy elements up to barium, and through the third ionization 
states. The line list was updated to include new wavelengths and/or oscillator strengths for 
~180 Fe III lines (Proflitt et al. 1999; Ekberg 1993). Brooks et al. also included some fine 



3 LINFOR was originally developed by H. Holweger, W. Steffen, and W. Steenbock at Kiel University. It 
has been upgraded and maintained by M. Lemke, with additional modifications by N. Przybilla 
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adjustments to iron-group features after an examination of the synthesis of their sharpest 
lined stars; this included 5 Mn III lines and 26 Fe III lines. The final linelist includes 6685 
features, although many are negligible contributors. Atomic data for the B III 2s 2 S - 2p 2 P 
resonance doublet with lines at 2065.8 A and 2067.3 A are taken from Profntt et al. (1999). 

The weaker B III line at 2067.3 A is blended with a strong Fe III line and weaker Mn III 
line, and is not suitable for boron abundance determinations. For all syntheses, an isotopic 
ratio of 11 B/ 10 B=4.0 is adopted, the solar system ratio (Zhai & Sahw 1994, Shima 1963) 
which is consistent with the values determined by Profntt et al. (1999) from their line profile 
analysis of two sharp-lined B-type stars. 

As discussed by Profntt & Quigley (2001), interstellar (IS) lines must be addressed, 
particularly due to the capacity of one Cr II IS line at A2065.50 to be shifted into a blend 
with the boron feature. All of our stars have been examined for IS lines. In no instance does 
an IS line interfere with the B III features. All four of the IS lines near the B III feature are 
only observable in one star, HD 22951 (see Fig. 4). 

3.1. Spectrum Syntheses: I. Line by Line Analysis 

For five of the seven targeted stars, the S/N of the reduced spectra and their broadening 
parameters (sharp lines) made it possible to identify and fit individual stellar features for 
a line-by-line analysis. The wavelengths and abundances from these (mostly iron-group) 
features are listed in Table 4. Abundance results are listed relative to meteoritic values 
from Asplund, Grevesse & Sauval (2005), with 12+log(Fe/H, Mn/H, Ni/H, Cr,H, Co/H) 
= (7.45, 5.47, 6.19, 5.64, 4.82). Not all iron-group features were used in the final iron-group 
abundances listed in Table 5; due to poor atomic data in the UV and/or unrecognized blends, 
we iterated such that the average value includes only features within 2a of the mean. 

Microturbulence values were calculated in this line-by-line analysis by requiring that 
iron-group abundances match for both strong and weak features (i.e., the features in Table 4 
within 2a of the mean). This analysis yields microturblence results significantly lower than 
ones determined using optical observations of these same stars (Gies & Lambert; 1992). 
This same trend has been noted by both Venn et al. (1996) and Cunha et al. (1997), and is 
likely a reflection of differences in the atmospheric structure between the formation of the 
UV and optical spectral regions. Macroturbulence was taken as the instrumental broadening 
(9 km s _1 for the E230M grating and 4 km s _1 for the E230H grating). Values for rotational 
broadening and radial velocities were calculated via line profile fitting for each star from 
the iron group features in Table 4. Spectrum syntheses in Fig. 4 show the best fit iron 
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abundance. 

Boron abundances were determined from best fit syntheses of the A2065.8 B III line, 
even though it is blended with a weak Mn III line at A2065.9. The measurement of this 
Mn III line has been shown to be in good agreement with the overall iron-group abundance 
(Venn et al. 2002 = V02), thus the mean Fe group abundance was adopted for the Mn 
abundance prior to the boron abundance calculations. The best fit boron abundances, and 
syntheses with [B/H]=±0.4 dex, are shown in Fig. 4. The broadening parameters, and LTE 
iron-group and boron abundances are summarized in Table 5. 

In this analyses, all of the target stars result in lower iron-group abundances than an- 
ticipated for objects in the Orion star forming region or the solar neighborhood. Cunha, 
Smith & Lambert (1998) examined F and G-type stars in the Orion complex and found iron 
abundances from optical observations of Fe I and Fe II lines that were consistent with other 
stars in the solar neighborhood (log(Fe/H)~ —0.1). The fact that our analysis of B-stars 
results in iron-group abundances that are sub-solar by [Fe/H] = 0.1 to 0.4 appears to be 
partially due to the temperatures adopted. In V02, boron and iron-group abundances were 
determined using temperatures from the literature. The boron abundances were then scaled 
to a common temperature scale, but not the iron-group abundances. If the iron-group abun- 
dances had also been scaled (e.g., reducing the Gies & Lambert 1992 temperatures by 3.4%, 
as done in this paper from the beginning), then those abundances would have been reduced 
by ~0.1 dex, in better agreement with the results in this paper. The iron-group abundances 
are also very sensitive to microturbulence for which there are very few constraints in our UV 
spectra. In V02, we showed that A£ = =(=1 km s _1 results in A [Fe/H] ~ ±0.20. If all of 
the microturbulence values in this paper are reduced by 1 km s~ 1 (with no changes in tem- 
perature) then the mean iron-group abundance is solar. These changes in temperature and 
microturbulence have much smaller effects on boron, A[B/H] <0.05 and 0.02, respectively 
(see Section 4). Thus, we do not expect that the low iron-group abundances found here 
affect the main science in this paper, i.e., accurate boron abundances in B-type stars. 



3.2. Spectrum Synthesis: II. Two Exceptions 

Of the seven target stars, two had spectra that prevented a detailed line-by-line study. 
Because only three of the subexposures for HD 214993 could be used, this resulted in a lower 
signal-to-noise spectrum, while HD 30836 has very broad features due to its high vsini value. 

Lines that were determined to be "clean" in the high quality spectra (lines that do 
not deviate from the mean abundance by more that 2a for the 5 stars discussed above in 
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Section 3.1, see Table 4) were used to find the rotational and radial velocities. Using these 
parameters, a best metallicity and microturbulence was determined simultaneously (where 
the metallicity range examined ranged from solar and 1/3 solar). The best-fit spectrum 
syntheses are shown in Fig. 5. In both cases, sub-solar metallicities were found to fit best 
(see Table 5). For HD 30836 and HD 214993 the best fit iron-group abundance was adopted 
for the Mn III abundance (blended with the boron line), and used to calculate the boron 
abundance or its upper limit (see Fig. 5). 



3.3. Spectrum Synthesis: III. HD35299 

We have synthesised the Goddard High Resolution Spectrograph (GHRS) spectrum of 
HD 35299 near 2066A from Proffitt et al. (1999; see also Lemke et al. 2000) to verify that 
our methods and models are consistent with previous analyses that found a solar boron 
abundance for this star. Temperature and gravity values are from CL94. Macroturbulence 
was set to the instrumental broadening, 3.5 km s -1 , as calculated using FWHM (0.0246 A) 
from Proffitt et a/.(1999). Rotational broadening (8 km s -1 ), and radial velocities were 
calculated here via line profile fitting of the iron features deemed best in Table 4 (although 
the spectral range is smaller AA2059-2073). Microturbulence was taken as the best fit value 
for each of the input abundances. The HD 35299 spectrum is plotted in Fig. 6 and shows the 
best-fit iron-group (log(Fe/H)=7.20) and boron (log(B/H)jvz,T£ = 2.55) abundance synthetic 
spectrum. For comparison, we also show the best-fit boron abundance ±0.4 dex. This 
iron-group abundance is consistent with the average of log(Fe/H)from the two "clean" (as 
described in Section 3.1) iron lines that fall into the wavelength range of our observations. 
This value was used to set the abundance of the Mn III blend at A2065.9. This NLTE boron 
abundance is in good agreement with the results from both Proffitt et al. (log(B/H)=2.41) 4 
and Lemke et al. (log(B/H)=2.70). 



3.4. Nitrogen and Oxygen Abundances 

In Table 3, we have reduced the GL92 temperatures by 3.4% and applied their nitrogen 
and oxygen abundance corrections for this temperature change (the A values in their Ta- 
ble 9). We have also applied a correction to account for differences in the LTE abundances 
determined from their use of Gold as opposed to the more heavily line-blanketed Kurucz 



4 This value has been adjusted to the CL94 temperature scale and NLTE corrected using Table 8 from 
V02 
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model atmospheres. These corrections have been calculated by CL94 (their Table 10). Note 
that CL94's corrections are to the LTE abundances, and we assume that the same will apply 
to the NLTE abundances. For stars that have been studied only by CL94, similar correc- 
tions have been applied to account for the differences between LTE Gold and Kurucz model 
abundances. 



4. Boron and Iron-group Abundance Uncertainties 

In this analysis, we have listed the line-to-line scatter in the iron-group abundances 
in Table 4. For all other iron-group uncertainties (atmospheric parameters, metallicities, 
spectral resolution, and S/N), we refer to Table 7 in V02 since their targets stars, UV 
spectra, and methodology are similar to those in this paper. Most of these atmospheric 
effects result in iron-group uncertainties <0.1 dex, with the exception of £ which is larger 
(~0.20 dex). The macroturbulence and radial velocity measurements are estimated with 
uncertainties of ±2 km s _1 based on the line profile fitting. Microturbulence is estimated as 
accurate to ± 1 km s' 1 based on line strengths. Overall, the iron-group uncertainties are 
estimated as Alog(Fe/H)=0.25 dex. 

To determine the uncertainties in the boron abundances, we have used similar meth- 
ods to those described in V02 and Brooks et al. 2002, calculating the offset in the boron 
abundance due to the uncertainties in temperature and gravity (as published), macroturbu- 
lence, microturbulence, and radial velocities (from our iron-group analysis), and location of 
the continuum (±1%). The abundance of the Mn III line blended with the B III line was 
maintained constant during the boron uncertainty calculations. Brooks et al. showed that 
the error determined when Mn remains fixed is in good agreement with that determined 
by refitting the Mn line with each successive error calculation. The resulting boron uncer- 
tainties are listed in Table 6; the total error for boron is taken as the sum of the squares 
of these values. As was found in Brooks et al, the accuracy of the atmospheric parameters 
(temperature and gravity) is the greatest source of error in the boron abundance. 

The isotopic ratio of 11 B/ 10 B can also have a small effect on the derived boron abun- 
dance. In V02, a ratio of 2.0 was found to change log(B/H) by <0.12 dex. 

4.1. The Lyubimkov et al. Temperature Scale 

In addition to the uncertainties in the atmospheric parameters, we also calculated the 
errors if we adopt a different temperature scale. Lyubimkov et al. (2002) report cooler 
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temperatures for B-stars. The iron-group and boron abundances were calculated for five of 
our stars (those that underwent the line-by-line analyses) with temperatures and gravities 
from both Lyubimkov et al. and GL92-3.4%, and are compared in Table 7. Both the iron 
and boron abundances are lower by 0.02 to 0.26 dex using the Lyubimkov temperatures (and 
gravities). There is no obvious correlation between these corrections and temperature, unlike 
the results for oxygen from Daflon et al. (2004). The typical correction to the iron or boron 
abundances is ~ —0.1 dex, thus the change in temperature scale seems to have only a small 
effect on our abundance determinations. 

We also note that for this test, four (out of the five) stars had values of T e g and 
log g calculated directly by Lyubimkov et al. (2002). In order to correct the fifth star, 
HD 22951, equations 6 and 8 from Lyubimkov et al. were used to calculate the new T e g. 
This choice of equations was based on a previous log g=4.4 determination for HD 22951. 
This new temperature was then used to estimate a new value for log g (from their Fig. 13). 
Values for [ci], Q, and /3 were obtained from Hauck & Mermilliod (1998). In order to confirm 
our methods, adjusted values of T c g and log g were calculated for HD 35337 with the above 
described method and were found to be in excellent agreement with the values listed in 
Table 1 of Lyubimkov et al. (2002). We note this because the iron and boron abundance 
corrections for HD 22951 are the largest in Table 7. 



5. Discussion 

Current stellar evolution models that include rotation follow the evolution of the angular 
momentum distribution and associated mixing processes in massive stars from the pre-main 
sequence through core collapse. Rotational mixing can affect stellar surface abundances, 
stellar life times, and the evolution of a star across the HR diagram, thus stellar mass 
estimates. Generally, a rotating star has a lower effective gravity, thus it acts like it has less 
mass at core-H ignition. During core-H burning, rotationally induced mixing of protons from 
the envelope into the convective core and of helium from the core into the envelope will lead 
to higher luminosities compared to non-rotating models. Evolution on the main-sequence 
then depends on mass and rotation rate, as well as the efficiency of mixing in the upper 
convective core and stellar interior. 

Stellar evolution models by Heger & Langer (2000) follow the changes in the surface 
abundances of boron and nitrogen as a function of time (age) for a variety of stellar masses 
(8 to 25 M Q ) and initial rotation rates (0 to ~450 km s" 1 ). These relationships are impor- 
tant because they show that boron is a much more sensitive indicator of mixing processes 
in B-type stars than the more commonly examined CNO abundances. Two sets of models 
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are computed by Heger & Langer, with different assumptions on the efficiency of rotational 
mixing in layers containing a gradient in the mean molecular weight \x (called a /i-barrier). 
The differences between the two sets of models reflects the remaining uncertainties in the 
theoretical description of the rotational mixing processes. It is interesting to note that even 
though the rotationally induced ^-barrier inhibits mixing just above the core, the current 
models show that the envelope above the ^-barrier is very well mixed, thus CN-processed gas 
is more evident up through the photosphere than when the /z-barrier is ignored. The models 
also show the degeneracies in the stellar abundances (CNO and boron) with mass, age, and 
rotation rate. Thus, determination of these stellar abundances alone cannot uniquely con- 
strain rotational mixing scenarios; additional information (such as mass, rotational velocity, 
or age) of the target stars is necessary. 

In this paper, we attempt to constrain both mass and age for our selection of stars while 
interpreting their CNO and boron abundances. Mass is constrained by the model atmosphere 
parameters; in Fig. 7, it is clear that most of our target stars are ~12 M , ranging in mass 
from 8 to 14 M Q . Age is constrained by our initial selection of stars in young clusters 5 , 
see Table 1. Unfortunately, the rotation rates themselves cannot be constrained because 
we are required to analyse sharp-lined stars for reliable boron determinations, which could 
be either slowly-rotating or nearly pole-on rapidly rotating stars. Abt, Levato, & Grosso's 
(2002) statistical survey of rotation rates of B-type stars in the Bright Star Catalogue found 
B0-2V stars with a mean vsini value of 127 ±8 km s" 1 , and that the highest vsini value is 410 
km s -1 , though these are very rare stars (0.3% of the sample had vsini > 350km s _1 ). This 
velocity range is covered by the Heger & Langer models. Nevertheless, with nearly all other 
parameters constrained, then we can compare the predictions from the evolution models for 
boron versus nitrogen and mass and age, for a range of rotation rates for our stellar sample. 



5.1. Initial Abundances of Boron and Nitrogen 

The meteoritic abundance of boron is adopted as the initial abundance of boron in the B- 
type stars in this paper, 12+log(B/H)= 2.78±0.04 (Zhai & Shaw, 1994). This value is in good 
agreement with the new 3D solar photospheric boron abundance determined by Asplund, 
Grevesse, & Sauval (2005), as 12+log(B/H) = 2.70±0.2, with corrections for departures from 



°Thc cluster ages depend on isochrone fitting, which in turn depend on stellar evolution models. Standard 
stellar evolution models without rotation are usually adopted when cluster ages are determined. Rotational 
mixing could have an effect on the ages of clusters, but it is not clear what the effect on age would be. 
Isochrone fitting takes into account all main-sequence (turn off) stars, and these stars can be expected to 
display a variety of rotation rates. 
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LTE. Their solar nitrogen abundance is determined from 3D analyses of NI and NH features 
(the NI result includes a small correction for NLTE), resulting in 12+log(N/H) = 7.78 ±0.06, 
which is in good agreement with other stars in the solar neighborhood. Additionally, the 
interstellar boron abundance from B II 1362 A resonance lines in diffuse clouds in the solar 
neighborhood has a lower limit of 12+log(B/H)> 2.4 (Howk, Sembach, & Savage 2000, where 
a lower limit is given since boron may be depleted onto interstellar grains). In this paper, we 
assume the solar abundances are representative of the stars in the solar neighborhood, even 
the younger B-type stars (even though one may ask why there has been so little chemical 
enrichment over the ~5 Gyr lifetime of the Sun). 



5.2. Boron versus Nitrogen: Predictions and Observations 

In hot stars, the boron-nitrogen relationship is of particular interest because the ma- 
jority of boron in the star is destroyed very early in its main-sequence lifetime (< 10 4 yr), 
leaving only the outermost layers to retain their initial boron abundance. Subsequent mixing 
of CN-cycled gas towards the surface of the star through convective overshoot and/or semi- 
convection results in a downward mixing of the boron rich outer layers toward the stellar 
interior. Because of the shallow mixing required for boron to be destroyed (only about 1 M Q 
down), its depletion manifests itself prior to the CN-cycle nitrogen gas being mixed to the 
surface. Predictions for the change in boron versus nitrogen from the Heger & Langer (2000) 
models are shown in Fig. 8, where it can be seen that a clear signature of rotating models is 
that the reduction in boron is far greater than the increase in nitrogen. The dotted lines in 
Fig. 8 are from the same model (12 M , ZAMS V rot = 200km s _1 , and including yu-barriers) , 
but various initial boron and nitrogen abundances. Interestingly, the boron-to-nitrogen re- 
lationship shown in Fig. 8 is nearly independent of mass, rotation velocity, and age in the 
stellar models (see Fig. 12 in V02). 

Fig. 8 also shows that that most of the observational data follows or exceeds the predicted 
trend. All filled data points in Fig. 8 are boron determinations from B III 2066. Only one 
star has a solar boron abundance as determined from the B III lines, HD 35299, thus we 
reanalysed it here (see Section 3.3) to verify our methods. All the other stars have boron 
results from the B III line below 12+log(B/H)=2.3. For the stars in this paper, lower boron 
abundances are to be expected since targets were selected to have weak B III lines from IUE 
spectra. 

Those stars with larger boron depletions than predicted from models may indicate addi- 
tional or more efficient mixing, but one of the major advantages of the data set in this paper 
is that the boron and nitrogen abundances have been determined for stars in clusters, thus 
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stars which also have age constraints. In Figs. 9 and 10, we show the nitrogen and boron 
abundances versus their stellar cluster ages, respectively. Ages are listed in Table 1. The 
models shown in Figs. 9 and 10 (and also Fig. 7) include rotational mixing and /z-barrier 
effects for several different stellar models across a span of rotational velocities (0-450 km s _1 ) 
from 10 to 20 M in order to illustrate the effects of these parameters on abundance. 

The predicted trend of boron-depletion and nitrogen-enrichment lends itself to certain 
constraints in the ability to detect boron if nitrogen is enhanced and/or with age. For 
example, boron is nearly undetectable after a nitrogen enrichment of only 0.4 dex (based 
on Fig. 8). From Fig.s 9 and 10, for a 12 M & star with v rot =400 km s -1 , this time scale is 
no more than 10 Myr. Thus, boron's sensitivity to early mixing is far greater than that of 
nitrogen. 

In comparing the observational stellar abundances to the model predictions, the stellar 
data do show a large range in boron abundances, from 0.9 > 12+log(B/H) > 2.9. The 
stars with 12+log(B/H) > 1.8 (which comprises 3/4 of the sample) do not show significant 
nitrogen enrichments; see Figs. 8 and 9. This suggests a range of a factor of 10 in boron 
without nitrogen enrichments, as predicted by stellar evolution models that include rotation. 
The remaining 1/4 of the sample, with 12+log(B/H) < 1.8, show a large range in nitrogen, 
some with no enrichments at all, i.e., HD 36591 and HD 30836. These two stars are discussed, 
independently, below. 

In Fig. 10 most of the stellar data is in good agreement with the model predictions (for 
~12 M stars, and given the uncertainty in the rotation rates). Only two stars stand out, 
HD 36591 and HD 205021, because their ages are so young and boron abundances so low. 
For HD 36591, its young age and low boron abundance suggests a 20 M model with the 
fastest rotation rate (450 km s _1 ); see Fig. 10. On the other hand, examination of Fig. 9 
suggests a 15 M Q model with the fastest rotation rate. However, its spectroscopic mass 
(from Fig. 7) is at most 13 M & . In order for the masses to agree, the temperature needs to 
be increased by ~15-30%. This result suggests that rotational mixing may be more efficient 
than currently modelled, at least at the highest rotation rates, and that the current models 
do not simultaneously predict the boron abundance, nitrogen abundance, mass, and age for 
this star. 

The remaining two stars in our sample that stand out in either Fig. 8 or Fig. 10 are 
HD 30836 and HD 205021. Both are severely depleted in boron, but have also been identified 
as spectroscopic binaries. Binary mass transfer is another way of depleting surface boron 
abundances (Wellstein 2001; Langer, Wellstein, & Heger 2001). These stars have similar 
masses (1O-12M from Fig. 7) and orbital periods (10 days, Koch 1990). Models from Well- 
stein (2001) predict that mass transfer occurs in these systems shortly after core hydrogen 
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burning in the primary star. While binary mass transfer can produce boron depletions, only 
rotational mixing predicts boron depletion without nitrogen enrichments. HD 205021 shows 
a nitrogen enrichment (log(N/H) =8.0) in addition to its depleted boron abundance. We 
notice in Fig. 9 that the nitrogen enrichment and age for this star implies a mass of ~20 
M Q and the highest rotation rate (450 km s _1 ). The same result is found examining its 
low boron abundance and age in Fig. 10. However, its spectroscopic mass is only ~12 M 
(Fig. 7). Thus this star does not match the rotating model predictions, but could be affected 
by binary mass transfer. 

HD 30836 shows no nitrogen enrichment (log(N/H)=7.79; Table 3), thus rotational mix- 
ing uniquely describes its abundance pattern. Ironically, we only have an upper-limit to the 
boron abundance in HD 30836, however that upper-limit is extremely low. We also only have 
an upper-limit to its age from membership in the Orion OBI association (the sub-association 
is unclear for this star). From the rotating models, the upper-limit age and upper-limit boron 
abundance suggest that HD 30836 has a mass of 12 M Q with a rotational velocity between 
200 and 300 km s -1 (see Fig. 10). However, the spectroscopic mass (see Fig. 7) for this star 
is closer to 10 M . Adopting 10 M in Fig. 10 would require a rotational velocity over 
450 km s _1 , higher than the theoretical upper limit. Similarly, adopting a 12 M & in Fig. 7 
would require an increase of nearly 3000 K (approximatly 15%) in temperature. Thus, the 
nitrogen abundance and upper-limits to age and boron abundance make HD 30836 another 
interesting star for comparing with the rotating stellar models, regardless of its binarity. 
This star, along with HD 36951 (and possibly HD 205021) are best fit by models with masses 
significantly larger than their spectroscopic masses. This suggests that rotational mixing is 
more efficient than currently modelled at the highest rotation rates. 



6. Conclusions 

Detailed abundances of boron have been determined in a careful selection of seven sharp- 
lined B-stars in young clusters to test massive star evolution scenarios that include rotational 
mixing effects. The seven stars in this analysis show moderate-to-severe depletions of boron 
(1/4 to 1/50 solar!) without significant nitrogen enrichments. This is a unique prediction 
of the massive star evolution models that include rotational mixing (Heger & Langer 2000), 
because shallow mixing will destroy surface boron abundances before deep mixing bring CN- 
cycled gas to the surface. Only three stars deviate from the model predictions, HD 30836, 
HD 36591, and HD 205021. In all three cases, their spectroscopic masses are smaller than 
predicted from the rotating evolution models (i.e., the masses required to explain their low 
boron abundances and young ages). These three stars appear to indicate that rotational 
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mixing is more efficient than currently modelled at the highest rotation rates. These results 
are consistent with previously published results (Venn et al. 2002), but approximately doubles 
the sample size of stars with both moderate and severe boron depletions. 

We would like to thank Katia Cunha for help in analysing HD 35299, and many pro- 
ductive and enjoyable conversations. This research was supported by NASA grant HST 
GO-09437. Space Telescope Science Institute is operated by the Association of Universi- 
ties for Research in Astronomy, Inc. under NASA contract NAS 5-26555. DLL thanks the 
Robert A. Welch Foundation of Houston, Texas for their support. 
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Table 1. Basic Parameters 



Assoc. Age 

Star Name Sp. Type Variable Binary OB Assoc. (Myr) Ref. 



HD 22951 


o Per 


B0.5V 


Visual? (check) 


Per 0B2 


12-20 


1,3,4 


HD 30836 


7T 4 Ori 


B2III 


Spect. 


Ori OBI 


<11.4 


2,5,6 


HD 34816 


A Lep 


B0.5IV 




Ori OBI 


<11.4 


2 


HD 35337 


8 Lep 


B2IV 




Ori OBlc 




2,6 


HD 36285 




B2IV-V 




Ori OBlc 




2,6 


HD 36960 




B0.5V 


Variable? Visual 


Ori OBlc 




2,7,8 


HD 214993 


12 Lac 


B2III 


Cep 


Lac OBI 


6-8 


1,9 



References. - - (l)de Zeeuw & Brand 1985 (2)Brown, de Geus, & de Zeeuw 1994 (3)Klockova & 
Kopylov 1985 (4)Petrova & Orlov 1999 (5)Batten, Fletcher, & Mann 1978 (6)Abt & Cardona 1984 
(7)Morrell & Levato 1991 (8)Lindroos 1986 (9)Struve 1951 
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Table 2. HST STIS Observing Information for Galactic B-stars 



Star 


V 


Grat/Slit 


Exposurc(s) 


Date 


S/N 


HD 22951 


4.98 


E230M 


1707 at A c 2124 


2003 Mar 25 


50 






0.2x0.05ND 


2634 at A c 2124 












2773 at A c 1978 












2760 at A c 2269 






HD 30836 


3.67 


E230M 


856 at A c 2124 


2003 Apr 6 


65 






0.2x0.05ND 


720 at A c 2124 












1200 at A c 1978 












1091 at A c 2269 












720 at A c 2124 


2003 Aug 2 










600 at A c 2124 












1200 at A c 1978 












1200 at A c 2269 






HD 34816 


4.27 


E230M 


855 at A c 2124 


2003 Apr 25 


80 






0.2x0.05ND 


720 at A c 2124 












1200 at A c 1978 












1091 at A c 2269 






HD 35337 


5.22 


E230M 


3036 at A c 2124 


2003 Apr 24 


100 






0.2x0.05ND 


2166 at A c 1978 












1347 at A c 2269 






HD 36285 


6.32 


E230H 


572 at A c 2013 


2003 Apr 28 


280 






31x.0.05NDA 


1080 at A c 2013 












1260 at A c 2013 












1 979 at \ 901 "} 












1200 at A c 2013 












624 at A c 2013 












1410 at A c 2013 












1410 at A c 2013 












1410 at A c 2013 












1388 at A c 2013 






HD 36960 


4.78 


E230M 


1661 at A c 2124 


2003 Apr 27 


90 






0.2x0.05ND 


1560 at A c 1978 












720 at A c 2269 






HD 214993 


5.23 


E230M 


2160 at A c 2124 


2003 Jul 7 


40 






0.2x0.05ND 


2640 at A c 1978 


2003 Jul 8 










1800 at A c 2269 
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Table 3. Stellar Parameters from Literature 





GL92-3.4% 


CL94 




NLTE 


NLTE 


Corrected 


Corrected 


Star 


T c tt 


T e s 


log g 


log(N/H) 


log(0/H) 


log(N/H) 


log(0/H) 




(K) 


(K) 












HD 22951 


27870 




4.40 


7.69 


8.45 


7.69 


8.42 


HD 30836 


21370 




3.60 


7.79 




7.75 




HD 34816 


28870 




4.20 


7.66 


8.75 


7.59 


8.67 


HD 35337 


23590 




4.20 


7.65 


8.56 


7.64 


8.55 


HD 36285 




21930 


4.40 


7.77 


8.79 


7.77 


8.55 


HD 36960 


28940 


28920 


4.30 


7.72 


8.88 


7.65 


8.80 


HD 214993 


24760 




4.00 


7.80 


8.78 


7.82 


8.83 



Note. — Corrections made to the nitrogen and oxygen abundances are based on both the reduced 
GL92 temperature scale and an adjustment from the less heavily line-blanketed Gold models to 
those of Kurucz, as discussed by CL94. Adjustments are made according to the errors found in 
Tabic 9 of GL92 (A values) and Table 10 of CL94. 
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Table 4. Iron-group Abundance Results 



A 




[M/H] 


[M/H] 


[M/H] 


[M/H] 


M/H] 


(A) 


Element (s) 


HD 22951 


HD 34816 


HD 35337 


HD 36285 


HD 36960 



2048.92 Mn III +0.26 

2049.37 Felll+MnUI -0.20 

2049.66 Mn III -0.51 

2050.74 FeUI -0.90 

2051.85 FcIII+FelV -0.05 

2052.27... FeUI -0.60 

2053.52 FeUI -1.00 

2054.56... Fe 111x3 -0.37 

2055.86 Fe III -0.17 

2056.16 Fe III -0.59 

2057.07 Fe III -0.54 

2057.93... Fe III -0.49 

2058.21 FeUI -0.11 

2058.57 Fe III -0.35 

2059.67 FeUI -0.82 

2063.40 Mn III -0.42 

2066.40... Mn IIIx2+Ni III -0.54 

2068.26 Fe III -0.57 

2068.99 Fe III+Mn Ill+Cr III -0.70 

2069.82... Fe III+Mn III -0.20 

2070.56 Fe III -0.40 

2070.98 Fe 111x3 -0.55 

2073.35... Fe III+Mn III -0.30 

2074.23 Fe III -0.62 

2076.32 Fe III -0.63 

2077.36... Mnlll+CoIII -0.23 

2077.74 Felll+FcIV +0.12 

2078.08 Fe III+Mn III -0.30 

2079.00 Fe 111x4 -0.60 

2080.22... FeUI -0.30 

2081.08... Mn IIIx2+Co III -0.25 

2082.38... FeUI -0.25 

2083.55... FeUI -0.20 

2084.36 Fe IIIx3+Mn 111x2 -0.60 

2084.93 Fe 111x2 -0.40 

2085.84... Felll+Crlll -0.25 



-0.09 


+0.13 


-0.04 


-0.06 


+0.00 


-0.08 


-0.05 


-0.10 


-0.43 


+0.38 


+0.16 


-0.31 


-0.25 


+0.01 


-0.35 


-0.12 


-0.15 


-0.09 


-0.10 


-0.15 


-0.41 


-0.17 


-0.35 


-0.27 


-0.71 


-0.43 


-0.38 


-0.45 


-0.23 


-0.13 


-0.04 


-0.10 


+0.02 


+0.03 


-0.21 


+0.14 


+0.13 


-0.20 


-0.25 


+0.15 


0.00 


-0.22 


-O.46 


+0.25 


-0.10 


-0.27 


-0.40 


-0.09 


+0.15 


-0.05 


-0.08 


+0.25 


+0.35 


+0.17 


-0.10 


+0.35 


+0.07 


+0.15 


-0.20 


-0.01 


-0.54 


+0.16 


-0.14 


-0.34 


-0.54 


-0.14 


-0.28 


-0.39 


+0.38 


+0.25 


+0.05 


+0.38 


-0.70 


-0.10 


-0.05 


-0.50 


-0.20 


-0.09 


-0.09 


-0.15 


-1-0 f)8 


-U. -L^; 


-u.ou 


-i-f) 1 ^ 


-0.48 


-0.32 


-0.32 


-0.33 


-0.28 


-0.05 


-0.13 


-0.30 


-0.65 


-0.32 


-0.55 


-0.40 


-0.74 


-0.05 


-0.50 


-0.72 


-0.53 


+0.04 


-0.13 


-0.33 


+0.10 


+0.15 


+0.01 


+0.10 


-0.40 


-O.43 


-0.50 


-0.43 


+0.25 


+0.09 


-0.45 


+0.20 


-0.38 


-0.10 


-0.25 


-0.29 


-0.13 


-0.06 


-0.28 


-0.23 


-0.35 


+0.00 


-0.10 


-0.36 


-0.20 


+0.05 


-0.25 


-0.07 


-0.43 


-0.15 


-0.65 


-0.37 


-0.30 


+0.15 


-0.20 


-0.15 


-0.15 


+0.02 


-0.14 


-0.15 
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Table 4 — Continued 



A 




[M/H] 


[M/H] 


[M/H] 


[M/H] 


M/H] 


(A) 


Element (s) 


HD 22951 


HD 34816 


HD 35337 


HD 36285 


HD 36960 


2087.15 


Fe III 


-1.05 


-0.31 


+0.10 


-0.37 


-0.38 


2087.93 


Fe III 


-0.80 


-0.51 


+0.13 


-0.40 


-0.44 


2089.12... 


Fe III 


-0.65 


-0.03 


+0.01 


-0.45 


+0.02 


2090.16 


Fe IIIx4+Mn 111x3 


-1.10 


-0.45 


+0.05 


-0.15 


-0.45 


2091.35 


Fe 111x2 


-0.65 


-0.46 


-0.32 


-0.48 


-0.40 


2092.97 


Fe III 


-0.20 


+0.00 


-0.05 


-0.22 


+0.15 


2093.51 


Fe III 


-0.25 


-0.58 


-0.08 


-0.25 


-0.13 


2095.66 


Fc 111x3 


-0.25 


-0.44 


-0.15 


-0.05 


-0.20 


2096.42 


Fe III 


-0.60 


-0.40 


-0.14 


-0.23 


-0.05 


2099.30 


Fc 111x2 


-0.50 


-0.55 


-0.15 


-0.45 


-0.42 


2101.04... 


Fe III+Mn III 


-0.50 


-0.30 


-0.07 


-0.18 


-0.20 


2103.74 


Fc 111x4 


-1.00 


-0.48 


-0.40 


-0.60 


-0.39 


2104.96 


Fe IIIx2+Cr III+Ni III 


-0.32 


-0.20 


+0.15 


+0.03 


-0.18 


2105.59 


Cr III 


-0.15 


-0.19 


-0.12 


-0.10 


-0.22 


2107.34 


Fe III 


-0.45 


-0.20 


-0.36 


-0.40 


+0.00 


2108.64... 


T ^ T T T - T T T 

Fe III+Mn III 


-0.42 


-0.43 


-0.05 


-0.07 


-0.30 


2111.80 


Fc 111x2 


-0.30 


-0.50 


-O.43 


-0.45 


-0.43 


2113.34 


Fe IIIx2+Mn III 


+0.00 


-0.50 


-0.04 


-0.30 


-0.18 


2113.83... 


Fc IIIx2+Cr III 


-0.40 


-0.35 


+0.02 


-0.23 


-0.10 


2114.34 


Fe IIIx2+Cr III 


-0.15 


-0.25 


+0.00 


-0.11 


-0.10 


2114.88 


Cr 111x2 


-0.10 


+0.00 


+0.06 


-0.05 


+0.08 


211 6 59 


Fc 111x2 


-0.95 


-0.86 


+0.07 


-O.48 


-0.72 


2117.55... 


Cr III 


-0.20 


-0.27 


-0.14 


-0.20 


-0.02 


2118.49... 


Fc IIIx2+Cr III 


-0.45 


-0.40 


-0.05 


-0.25 


-0.30 


2120.77... 


Fe IIIx4+Fc IV 


-0.53 


-0.20 


-0.07 


-0.07 


-0.15 


2123.59... 


Fe IIIx3+Cr III 


-0.50 


-0.40 


-0.02 


-0.05 


-0.27 


2124.16 


Fe IIIx5+Co III 


-0.80 


-0.70 


-O.4O 


-0.31 


-0.65 


2125.18 


Fe IIIx5+Fc IV+Mn III 


-0.55 


-0.45 


-0.39 


-0.33 


-0.35 


2126.14... 


Mn III 


-0.28 


-0.43 


+0.00 


-0.04 


-0.39 


2129.68... 


Fc 111x3 


-0.20 


-0.27 


-0.08 


-0.04 


-0.23 


2134.83 


Fe 111x2 


-0.08 


-0.04 


+0.29 


+0.03 


+0.12 


2136.36 


Fc 111x2 


-0.56 


-0.74 


-O.49 


-0.65 


-0.65 


Average.... 




-0.42 


-0.33 


-0.07 


-0.22 


-0.23 


1 a 




+0.15 


+0.16 


+0.15 


+0.14 


+0.15 
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Note. — Results that fall > 2a from the mean are italicized and not included in the final average or 1 a 
values. Wavelengths in bold are those that show consistent results across all five stars (none fall outside 
2a). 
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Table 5. Atmospheric Analysis Results 





vsmi 






LTE 


NLTE 1 


Star 


km s^ 1 


km s _1 


log(Fe/H) 


log(B/H) 


log(B/H) 


HD 22951 


23 


1 


7.03±0.10 


2.01±0.12 


2.01 


HD 30836 


43 


4 


7.25±0.25 


<1.2 


<1.0 


HD 34816 


35 


3 


7.12±0.11 


2.17±0.13 


2.17 


HD 35337 


15 


1 


7.38±0.10 


2.21±0.12 


2.11 


HD 36285 


10 


2 


7.23±0.09 


2.02±0.15 


1.82 


HD 36960 


33 


3 


7.22±0.10 


1.81±0.16 


1.81 


HD 214993 


39 


5 


7.25±0.25 


2.2±0.2 


2.1 



Note. — Macroturbulence values were set by the instrumental broaden- 
ing; 9 and 4 km s _1 for observations with the E230M and E230H gratings, 
respectively. 

: NLTE corrections were mading according to the grid used in V02 (their 
Fig. 9) 



Table 6. Boron Abundance Uncertainties 





HD22951 


HD34816 


HD35337 


HD36285 


HD36960 


Parameters 


Alog(B/H)) 


Alog(B/H) 


Alog(B/H) 


Alog(B/H) 


Alog(B/H) 


AT e// = ±750 K 


±0.07 


±0.07 


±0.06 


T0.05 


±0.12 


A log g = ±0.2 


±0.08 


±0.07 


±0.10 


±0.13 


±0.04 


A( = ± 1 km s _1 


T-0.01 


T0.03 


T0.02 


T0.02 


TO. 01 


A£ mQ = ± 2 km s" 1 


±0.00 


±0.02 


±0.00 


±0.00 


±0.00 


Shift continuum ±1% 


T0.05 


T0.08 


T0.03 


T0.04 


T0.09 


AV rci d = ± 2 km s" 1 


±0.00 


±0.01 


±0.02 


±0.02 


±0.01 


Total Error 1 


±0.12 


±0.13 


±0.12 


±0.15 


±0.16 



Total error is taken as the sum of squares of all of the uncertainties. 
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Table 7. Atmospheric Parameters from Lyubimkov et al. (2002) 





Lyubimkov02 


A 


A 


Star 


TeS log g 


[Fe/H] 


[B/H] 


HD 22951 


26300 4.00 


0.26 


0.18 


HD 34816 


27900 4.20 


0.10 


0.11 


HD 35337 


22300 4.00 


0.12 


0.06 


HD 36285 


21300 4.25 


0.02 


0.05 


HD 36960 


27500 4.15 


0.15 


0.13 



Note. — A = log(X/H) adopteti - log(X/H) L02 
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52828.00 52828.05 52828.10 52828.15 

Observation date (MJD) 

Fig. 1. — Measured radial velocity of each STIS subexposure of the /3-Cephei star 12 Lac 
(HD 214993) is plotted as a function of the obervation start time. We have overplotted a 
velocity curve that adopts the pulsational frequencies and amplitudes found Mathias et al. 
(1994), but have adjusted the phases of each of the six pulsation modes to optimize the fit. 
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2045 2050 2055 2060 2065 2070 2075 

Wavelength (A) 

Fig. 2. — We show the individual spectra from order 99 of the STIS E230M grating for 
each individual subexposure of 12 Lac (HD 214993). The first subexposure is plotted at the 
bottom of the figure, and each of the subsequent spectra is offset by 3 x 1CT 10 flux units. 
The spectra are aligned to the laboratory wavelength scale in air. 
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Coaddition of subexposures 1—3 & 7—11 




Coaddition of subexposures 4, 5, & 6 



2045 2050 2055 2060 2065 2070 

Wavelength (A) 



2075 



Fig. 3. — Coadded spectra for the three most narrow-lined observations of 12 Lac (bottom 
panel) are compared to the coaddition of the data for other phases. Each subexposure was 
shifted to the laboratory wavelength scale before coadding. Wavelengths are given in air. 
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Fig. 4. — Best-fit metallicites and boron abundances (red) are compared to the spectra of five 
stars (black) with sufficiently high signal-to-noise spectra for a line-by-line analysis. Best-fit 
boron abundances ±0.40 dex are shown (blue) for comparison. Four interstellar lines are 
marked above the figures, while the boron lines are indicated below. 
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Fig. 5. — Best-fit metallicities and boron abundances (red) are compared to the spectra 
of two stars (black) with either lower signal-to-noise spectra (HD 214993) or lower gravity 
(HD 30836) than the other stars in this program. Best-fit boron abundances ±0.40 dex 
are shown (blue) for comparison, which also emphasizes the upper-limit determination for 
HD 30836. Rotational and radial velocities were determined from the most reliable iron- 
group lines (in Table 4 and marked above the spectra). Boron is indicated below the plots 
by the dotted lines. 
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Fig. 6.— The HST-GHRS spectrum of HD 35299 (black; Proffitt et al. 1999) with best fit 
iron-group (log(Fe/H)=7.20) and boron (log(B/H)=2.55) spectrum synthesis (red). Best-fit 
boron ±0.40 dex is also shown (blue). We use this star to check our methods (see Section 
3.3). 
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Fig. 7. — Temperature vs gravity to determine spectroscopic masses. Models with low (black 
line) to moderate rotation rates (dotted black line) and including /z-barriers (red line) are 
from Heger & Langer (2000). The effects of very rapid rotation on main-sequence position 
and thus spectroscopic mass can also be seen (green lines). Stellar masses associated with 
each model are labelled along the bottom of the figure. Only stars analysed in this paper 
and V02 are labelled. The temperatures of the stars have been scaled from their references 
to the Cunha & Lambert (1994) scale (see V02 for details). 
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Fig. 8. — Boron and nitrogen abundances in B-stars. Note that these abundances have 
been corrected from their original references to a common temperature scale (see V02 for 
details). One rotating stellar evolution model is shown from Heger & Langer (2000); that 
for 12 M , V ro t=200 km s _1 , and including the ^u-barrier (other masses and rotation rates 
have nearly identical predictions for the change in boron with nitrogen). The three dotted 
lines follow the predictions from this model for three different initial nitrogen and boron 
abundances; (12+log(N/H), 12+log(B/H)) = (7.6, 2.6), (7.7, 2.7), (7.8, 2.8). The x-axis and 
y-axis are scaled equally to emphasize the destruction of boron before significant nitrogen 
enrichment. 
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Fig. 9. — Nitrogen enrichment as a function of age. Data symbols the same as in Fig. 7. 
Rotating stellar evolution models are from Heger & Langer (2000). Red dashed lines rep- 
resent 12M Q models with various rotational velocities (0, 100, 200, 300, 450 km s _1 , where 
the lines for the slowest rotation rates lie along the bottom of the plot). Solid black lines 
represent 15 and 2OM models with the highest rotation rate (450 km s _1 ), and a 10M Q 
model with a moderate rotation rate (200 km s _1 ). 
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Fig. 10. — Boron depletion as a function of age. Data symbols the same as in Fig. 7. Rotating 
stellar evolution models are from Heger & Langer (2000). Red dashed lines represent 12M Q 
models with various rotational velocities (0, 100, 200, 300, 450 km s _1 ). Solid black lines 
represent 15 and 2OM models with the highest rotation rate (450 km s _1 ), and a 10M Q 
model with a moderate rotation rate (200 km s _1 ). 



